Effect of significant data loss on identifying electric signals that precede rupture by 

detrended fluctuation analysis in natural time 
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Electric field variations that appear before rupture have been recently studied by employing 
the detrended fluctuation analysis (DFA) as a scaling method to quantify long-range temporal 
correlations. These studies revealed that seismic electric signals (SES) activities exhibit a scale 
invariant feature with an exponent udfa ~ 1 over all scales investigated (around five orders of 
magnitude). Here, we study what happens upon significant data loss, which is a question of primary 
practical importance, and show that the DFA applied to the natural time representation of the 
remaining data still reveals for SES activities an exponent close to 1.0, which markedly exceeds the 
exponent found in artificial (man-made) noises. This, in combination with natural time analysis, 
enables the identification of a SES activity with probability 75% even after a significant (70%) data 
loss. The probability increases to 90% or larger for 50% data loss. 
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Complex systems usually exhibit scale- 
invariant features characterized by long-range 
power-law correlations, which are often diffi- 
cult to quantify due to various types of non- 
stationarities observed in the signals emitted. 
This also happens when monitoring geoelectric 
field changes aiming at detecting Seismic Electric 
Signals (SES) activities that appear before major 
earthquakes. To overcome this difficulty the 
novel method of detrended fluctuation analysis 
(DFA) has been employed, which when combined 
with a newly introduced time domain termed 
natural time, allows a reliable distinction of 
true SES activities from artificial (man-made) 
noises. This is so, because the SES activities 
are characterized by infinitely ranged temporal 
correlations (thus resulting in DFA exponents 
close to unity) while the artificial noises are not. 
The analysis of SES observations often meet the 
difficulty of significant data loss caused either by 
failure of the data collection system or by removal 
of seriously noise-contaminated data segments. 
Thus, here we focus on the study of the effect of 
significant data loss on the long-range correlated 
SES activities quantified by DFA. We find that 
the remaining data, even after a considerable 
percentage of data loss (which may reach ~ 80%), 
may be correctly interpreted, thus revealing 
the scaling properties of SES activities. This is 
achieved, by applying DFA not to the original 
time series of the remaining data but to those 
resulted when employing natural time. 



I. INTRODUCTION 

The output signals of complex systems exhibit fluctu- 
ations over multiple scalesii^ which are characterized by 



absence of dynamic scale, i.e., scale-invariant behavior—. 
These signals, due to the nonlinear mechanisms control- 
ling the underlying interactions, are also typically non- 
stationary and their reliable analysis cannot be achieved 
by traditional methods, e.g., power-spectrum and auto- 
correlation analysis^"—. On the other hand, the de- 
trended fluctuation analysis (DFA)ii^ has been estab- 
lished as a robust method suitable for detecting long- 
range power-law correlations embedded in non-stationary 
signals. This is so, because a power spectrum calculation 
assumes that the signal is stationary and hence when ap- 
plied to non-stationary time series it can lead to mislead- 
ing results. Thus, a power spectrum analysis should be 
necessarily preceded by a test for the stationarity of the 
(portions of the) data analyzed. As for the DFA, it can 
determine the (mono) fractal scaling properties (see be- 
low) even in non-stationary time series, and can avoid, in 
principle, spurious detection of correlations that are ar- 
tifacts of non-stationarities. DFA has been applied with 
successful results to diverse fields where scale-invariant 
behavior emerges, such as DNA^rJ^, heart dynamicsi^r— , 
circadian rhythms^^"— , meteorology^ and climate tem- 
perature fluctuations^"—, economics^"— as well as in 
the low- frequency (< 1 Hz) variations of the electric field 
of the earth that precede earthquakes^"— termed Seis- 
mic Electric Signals^^"^'* and in the relevant^"— mag- 
netic field variations^''. 

Monofractal signals are homogeneous in the sense that 
they have the same scaling properties, characterized lo- 
cally by a single singularity exponent ho , throughout the 
signal. Thus, monofractal signals can be indexed by a 
single global exponent, e.g., the Hurst exponent H = ho, 
which suggests that they are stationary from viewpoint 
of their local scaling properties (see Ref.^^- and refer- 
ences therein). Since DFA can measure only one expo- 
nent, this method is more suitable for the investigation 
of monofractal signals. In several cases, however, the 
records cannot be accounted for by a single scaling expo- 
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nent (i.e., do not exhibit a simple monofractal behavior). 
In some examples, there exist crossover (time-) scales 
separating regimes with different scaling exponents. In 
general, if a multitude of scaling exponents is required 
for a full description of the scaling behavior, a multi- 
fractal analysis must be applied. Multifractal signals 
are intrinsically more complex, and inhomogeneous, than 
monofractals (see Ref.^^ and references therein). A reli- 
able multifractal analysis can be performed by the Mul- 
tifractal Detrended Fluctuation Analysis, MF-DFA^i^ 
or by the wavelet transform (e.g., see Refj^). 

DFA has been applied, as mentioned, to the SES ac- 
tivities. It was found^ that when DFA is applied to 
the original time series of the SES activities and artifi- 
cial (man-made) noises, both types of signals lead to a 
slope at short times (i.e.. At < 30s) lying in the range 
a=l.l-1.4, while for longer times the range a^O. 8-1.0 
was determined without, however, any safe classification 
between SES activities and artificial noises. On the other 
hand, when employing natural time (see Section II) , DFA 
enables the distinction between SES activities and arti- 
ficial noises in view of the following difference: for the 
SES activities the a- values lie approximately in the range 
0.9 - 1.0 (or between 0.85 to 1.1, if a reasonable experi- 
mental error is envisaged), while for the artificial noises 
the a-values are markedly smaller, i.e., a=0.65-0.8. In 
addition, MF-DFA has been use d^^'^^ and it was found 
that this multifractal analysis, when carried out in the 
conventional time frame, did not lead to any distinction 
between these two types of signals, but does so, if the 
analysis is made in the natural time domain. 

The aforementioned findings of DFA for SES activi- 
ties, are consistent with their generation mechanising'^ 
which could be summarized as follows. Beyond the usual 
intrinsic lattice defects--"-- that exist in solids, in ionic 
solids in particular, when doped with aliovalent impuri- 
ties, extrinsic defects are formed for the sake of charge 
compensation. A portion of these defects are attracted 
by the nearby impurities, thus forming electric dipoles 
the orientation of which can change by means of a defect 
migration. In the focal area of an impending earthquake 
the stress gradually increases and hence affects the ther- 
modynamic parameters of this migration, thus it may re- 
sult in a gradual decrease of their relaxation time. When 
the stress (pressure) reaches a critical value, a coopera- 
tive orientation of these dipoles occurs, which leads to 
the emission of a transient signal. This signal constitutes 
the SES activity and, since it is characterized by critical 
dynamics, should exhibit infinitely range temporal cor- 
relations. This is consistent with the above findings of 
DFA that for SES activities udfa ~ 1- 

Hereafter, we will solely use DFA in view of its simplic- 
ity and its ability to reliably classify SES activities. It is 
the basic aim of this study to investigate how significant 
data loss affects the scaling behavior of long-range corre- 
lated SES activities inspired from the new segmentation 
approach introduced recently by Ma et al^ to generate 
surrogate signals by randomly removing data segments 



from stationary signals with different types of long-range 
correlations. The practical importance of this study be- 
comes very clear upon considering that such a data loss is 
inevitable mainly due to the following two reasons: First, 
failure of the measuring system in the field station, in- 
cluding the electric measuring dipoles, electronics and 
the data collection system, may occur especially due to 
lightning. Second, noise-contaminated data segments are 
often unavoidable due to natural changes such as rainfall, 
lightning, induction of geomagnetic field variations and 
ocean-earth tides besides the noise from artificial (man- 
made) sources including the leakage currents from DC 
driven trains. The latter are common in Japan where 
at some sites they may last for almost 70% of the time 
every day. We clarify, however, that even at such noisy- 
stations in Japan, several clear SES activities have been 
unambiguously identified^i during the night (when the 
noise level is low). In addition, prominent SES activi- 
ties were recently reported''^ at noise-free stations (far 
from industrialized regions) having long duration, i.e., of 
the order of several weeks. As we shall see, our results 
described in Section IV, are in essential agreement with 
those obtained in the innovative and exhaustive study 
of Ma et al^. Before proceeding to our results, we will 
briefly summarize DFA and natural time analysis in Sec- 
tion II, and then present in Section III the most recent 
SES data along with their analysis in natural time. In 
Section V, wc summarize our conclusions. 

II. CONVENTIONAL DETRENDED 
FLUCTUATION ANALYSIS. NATURAL TIME. 

We first sum up the original time series and deter- 
mine the profile y{i),i — 1, . . . ,N; We then divide this 
profile of length N into N/l{= Ni) non overlapping frag- 
ments of /-observations. Next, we define the detrended 
process yi^iy{m), in the ly-th fragment, as the difference 
between the original value of the profile and the local 
(linear) trend. We then calculate the mean variance of 
the detrended process: 

where 

/^(Z,^) = y ^y,%(m) (2) 

m— 1 

If F{1) P, the slope of the logF(/) versus log/ plot, 
leads to the value of the exponent a^pA = ol- (This 
scaling exponent is a self-similarity parameter that rep- 
resents the long-range power-law correlations of the sig- 
nal). If QfDi? A =0.5, there is no correlation and the signal 
is uncorrelated (white noise); if oldpa <0.5, the signal 
is anti-correlated; if oljjfa >0.5, the signal is correlated 
and specifically the case a£)i?^=1.5 corresponds to the 
Brownian motion (integrated white noise). 
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FIG. 1: (color online) Examples of the electric field recordings in normalized units, i.e., by subtracting the mean value fi and 
dividing by the standard deviation a. The following SES activities are depicted: (a)the one recorded on April 18, 1995 at 
loannina station, (b)the long duration SES activity recorded from December 27, 2010 to December 30, 2009 at Lamia station, 
(c) is an excerpt of (b) showing that, after long periods of quiescence, the electric field exhibits measurable excursions (transient 
pulses). 
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FIG. 2: (color online) (a) :Example of a surrogate time-series (in normalized units as in FiglTJ obtained by removing segments of 
length L = 200 from the signal of FiglTJa) with 50% data loss (i.e., p — 0.50). (b):The natural time representation of (a). The 
values obtained from the analysis of (b) in natural time are ki = 0.067(4), S = 0.076(4), S- — 0.071(4) and udfa ~ 0.90(5). 
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FIG. 3: (color online) The dependence of the DFA measure 
F{V] versus the scale / in natural time: we increase the per- 
centage of data loss p by removing segments of length L = 200 
samples from the signal of Fig[TJa). The black (plus) sym- 
bols correspond to no data loss (p=0), the red (crosses) to 
30% data loss (p=0.3), the green (asterisks) to 50% data loss 
(p=0.5) and the blue (squares) to 70% data loss (p=0.7). Ex- 
cept for the case p=0, the data have been shifted vertically for 
the sake of clarity. The slopes of the corresponding straight 
lines that fit the data lead to o-dfa =0.95, 0.94, 0.88 and 
0.84 from the top to bottom, respectively. They correspond 
to the average values of O-ofa obtained from 5000 surrogate 
time-series that were generated with the method of surrogate 
by Ma et al— (see the text). 
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We now summarize the background of natural time 
X- In a time series comprising TV events, the natural 
time Xk = k/N serves as an inde>s2^ for the occurrence 
of the k-th event. The evolution of the pair {xk,Qk) 
is studie d"''^ ,46/3- _78 ^ where Qk denotes a quantity pro- 
portional to the energy released in the fc-th event. For 
dichotomous signals, which is frequently the case of SES 
activities, the quantity Qk stands for the duration of the 
k-th pulse. By defining pk = Qk/Yln=iQn, we have 
found^'' that the variance 



(3) 



where (/(x)) = Y.k=iPkf{Xk), of the natural time x 
with respect to the distribution pk may be used^i^ for 
the identification of SES activities. In particular, the 
following relation should hold 



Ki « 0.070. 



(4) 



The entropy S in the natural time-domain is defined 

„46 



S" = (xlnx) - ix) ln(x>- 



(5) 



It exhibits^^ Lesch o^^i^° (experimental) stability, and for 
SES activities (critical dynamics) is smaller— than the 
value Su{= In 2/2 - 1/4 w 0.0966) of a "uniform" (u) 
distribution (as defined in Refs. — i^i^, e.g. when all pk 
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FIG. 4: (color online) The probabilities pi(a), P2(b) and P3(c) 
to recognize the signal of Fig[lja) as true SES activity when 
considering various percentages of data loss p=0.2, 0.3, 0.5, 
0.7 and 0.8 as a function of the length L of the contiguous 
samples removed. The removal of large segments leads to 
better results when using DFA in natural time (a), whereas 
the opposite holds when using the conditions of Eqs.(|6]) and 
([7]) for Ki, S and S- (b). The optimum selection (c) for the 
identification of a signal as SES activity consists of a proper 
combination of the aforementioned procedures in (a) and (b) , 
see the text. The values presented have been obtained from 
5000 surrogate time-series (for a given value of p and L) , and 
hence they a have plausible error 1.4% 1/^5000). 
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are equal or are positive independent and identically 
distributed random variables of finite variance. In this 
case, Ki and S are designated Ku{= 1/12) and Su, re- 
spectively.). Thus, S < Su- The same holds for the value 
of the entropy obtained '''^'''^ upon considering the time re- 
versal T, (the operator T is defined by Tpk = PN-k+i), 
which is labelled by S^. 

In summary, the SES activities, in contrast to the 
signals produced by man-made electrical sources, when 
analyzed in natural time (see Section I) exhibit in- 
finitely ranged temporal correlation al^' ^'^ and obey the 
conditions^!: 

Ki « 0.07 (6) 

and 

S,S-<Su. (7) 

It should be recalled that SES activities are publicized 
only when the magnitude of the impending earthquake is 
Ms(ATH)>6.0 

III. THE EXPERIMENTAL DATA 

In Figllja), we depict the SES activity recorded at 
loannina station (Northwestern Greece) on 18 April, 
1995. It preceded the 6.6 earthquake on 13 May, 1995. 
Since this earthquake was the strongest one in Greece 
during the 25 year period 1983-2007, we focus on this 
example in the next Section to present the DFA results. 

In addition, the most recent SES activity in Greece 
is depicted in FigslUb),(c). This has been recorded 
at Lamia station located in Central Greece during the 
period 27 December - 30 December, 2009. Almost 
three weeks later, two strong earthquakes of magnitude 
Ms{ATir)=5.7 and 5.6 occurred in central Greece with 
an epicenter at 38.4°N 22.0°E (but see also Refs,2ii2^). 

The two signals in Figs. [IIa),(b) have been classified 
as SES activities after analyzing them in natural time. 
In particular, for the signal in Fig. [IJa), straightfor- 
ward application of natural time analysis leads to the 
conclusion that the conditions (O and ([7]) are satisfied 
(see Table I of Ref;^). Further, the DFA analysis of the 
natural representation of this signal gives an exponent 
o^DFA = 0.95(4) « 1. (We also note that the classifica- 
tion of this signal as SES activity had been previously 
achieved by independent procedures discussed*"^ in the 
Royal Society Meeting that was held during May 11-12, 
1995 beforeS^ the occurrence of the 6.6 earthquake on 
May 13, 1995.) For the long duration signal of Fig. [^b), 
the procedure explained in detail in Ref4i was followed. 

IV. DATA ANALYSIS AND RESULTS 

Following Ma et al2S, we now describe the segmenta- 
tion approach used here to generate surrogate signals u{i) 



by randomly removing data segments of length L from 
the original signal u{i). The percentage p of the data loss, 
i.e., the percentage of the data removed, also character- 
izes the signal u{i). The procedure followed is based on 
the construction of a binary time-series g{i) of the same 
length as u(i). The values of u{i) that correspond to g{i) 
equal to unity are kept, whereas the data of u{i) when 
g{i) equals zero are removed. The values of u{i) kept, are 
then concatenated to construct u{i). 

The binary time-series g{i) is obtained as follows^S; 
(i)We first generate the lengths Ij — L with j = 
1, 2, . . . , M of the removed segments, by selecting M to 
be the smallest integer so that the total number of re- 
moved data satisfies the condition J2jLi h ^ ■ (ii) 
We then construct an auxiliary time-series a{k) with 
a{k) ~ L when k — 1,2,...,M and a(fe) = 1 when 
k ^ M + l,...,N - M{L + 1) of size N - M{L + 1). 
(iii)We shuffle the time-series a{k) randomly to obtain 
a{k). (iv)We then append a{k) to obtain g{i): if a{k) = 1 
we keep it, but we replace all a(k) — L with L elements 
of value '0' and one element with value '1'. In this way, 
a binary series g{i) is obtained, which has a size equal 
to the one of the original signal u{i). We then construct 
the surrogate signal u{i) by simultaneously scanning the 
original signal u{i) and the binary series g{i), removing 
the i-th element of u{i) if g(i) = and concatenating the 
segments of the remaining data to u(i). 

The resulting signal u{i) is later analyzed in natural 
time, thus leading to the quantities ki, S and S- as well 
as to the DFA exponent udfa in natural time. Such an 
example is given in Fig. [5J In what remains we present 
the results focusing hereafter, as mentioned, on the ex- 
ample of the SES activity depicted in Fig. [TJa). 

Typical DFA plots, obtained for L=200 and p = 30, 50 
and 70% are given in Fig. [31 For the sake of compari- 
son, this figure also includes the case of no data loss (i.e., 
p — Q). We notice a gradual decrease oi udfa upon in- 
creasing the data loss, which affects our ability to classify 
a signal as SES activity. 

In order to quantify, in general, our ability to identify 
SES activities from the natural time analysis of surrogate 
signals with various levels of data loss, three procedures 
have been attempted: 

Let us call procedure 1, the investigation whether 
otDFA^ resulted from the DFA analysis of the natural 
time representation of a signal, belongs to the range 
0.85 < aoFA < 1-0. If it does, the signal is then 
classified as SES activity. Figure SJa) shows that for a 
given amount of data loss (p=const), upon increasing the 
length L of the randomly removed segment, the proba- 
bility pi of achieving, after making 5000 attempts (for a 
given value of p and L), the identification of the signal 
as SES activity is found to gradually increase versus L 
at small scales and stabilizes at large scales. For exam- 
ple, when considering the case of 70% data loss (magenta 
color in Fig. Ufa)) the probability pi is close to 20% for 
L=50; it increases to pi ~ 30% for i=100 and finally 
stabilizes around 50% for lengths L =300 to 500. 
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Let us now label as procedure 2, the investigation 
whether the quantities ki, S and S'_ (resulted from the 
analysis of a signal in natural time) obey the conditions 
® and dZl), i.e., |ki-0.070| < 0.01 and S, S- <0.0966. If 
they do so, the signal is classified as SES activity. Figure 
IHJb) shows that for a given amount of data loss, the prob- 
ability p2 of achieving the signal identification as SES ac- 
tivity -that results after making 5000 attempts for each 
L value- gradually decreases when moving from the small 
to large scales. Note that for the smallest length scale in- 
vestigated, i.e., L=10 (which is more or less comparable 
-if we consider the sampling frequency of 1 sample / sec- 
with the average duration «11 sec of the transient pulses 
that constitute the signal) , the probability p2 reaches val- 
ues close to 100% even for the extreme data loss of 80%. 
This is understood in the context that the quantities ki, 
S and S- remain almost unaffected when randomly re- 
moving segments with lengths comparable to the aver- 
age pulse's duration. This is consistent with our earlier 
finding2^ that the quantities ki, S and S- are experimen- 
tally stable (Lesche's stability) meaning that they exhibit 
only slight variations when deleting (due to experimen- 
tal errors) a small number of pulses. On the other hand, 
at large scales of L, p2 markedly decreases. This may 
be understood if we consider that, at such scales, each 
segment of contiguous L samples removed, comprises on 
the average a considerable number of pulses the removal 
of which may seriously affect the quantities ki, S and 
As an example, for 80% data loss (cyan curve in 
Fig. IHb)), and for lengths L=400-500, the probability 
P2 of identifying a true SES activity is around 40%. 

Interestingly, a closer inspection of the two figures|4Ua) 
andlDJb) reveals that pi and p2 play complementary roles. 
In particular, at small scales of L, pi increases but p2 de- 
creases versus L. At large scales, where pi reaches (for 
considerable values of data loss) its largest value, the p2 
value becomes small. Inspired from this complementary 
behavior of pi and p2 , we proceeded to the investigation 
of a combined procedure, let us call it procedure 3. In 
this procedure, a signal is identified as SES activity when 
either the condition 0.85 < aoFA < 1-0 or the relations 
© and ([7]) are satisfied. The probability ps of achieving 
such an identification, after making 5000 attempts (for 
a given value of p and L), is plotted in Fig. Hl^c). The 
results are remarkable since, even at significant values 
of data loss, e.g., p = 70% or 80%, the probability p^ of 



identifying a SES activity at scales L=100 to 400 remains 
relatively high, i.e., ps w 75% and 65%, respectively (cf. 
note also that the value of reaches values close to 100% 
at small scales L=10). This is important from practical 
point of view, because it states for example the following: 
Even if the records of a station are contaminated by con- 
siderable noise, say 70% of the time of its operation, the 
remaining 30% of the non-contaminated segments have 
a chance of ~ 75% to correctly identify a SES activity. 
The chances increase considerably, i.e., to p^ ~ 90%, if 
only half of the recordings are noisy. 

The aforementioned results have been deduced from 
the analysis of a SES activity lasting around three hours. 
In cases of SES activities with appreciably longer dura- 
tion, e.g., a few to several days*^^!^ detected in Greece 
or a few months in Japan^^, the results should become 
appreciably better. 



V. CONCLUSIONS 

We start our conclusions by recalling that the distinc- 
tion between SES activities (critical dynamics, infinitely 
ranged temporal correlations) and artificial (man-made) 
noise remains an extremely difficult task, even without 
any data loss, when solely focusing on the original time 
series of electrical records which are, of course, in conven- 
tional time. On the other hand, when combining natural 
time with DFA analysis, such a distinction becomes pos- 
sible even after significant data loss. In particular we 
showed for example that even when randomly removing 
50% of the data, we have a probability (ps) around 90%, 
or larger, to identify correctly a SES activity. This prob- 
ability becomes somewhat smaller, i.e., 75%, when the 
data loss increases to 70%. To achieve this goal, the 
proper procedure is the following: the signal is first rep- 
resented in natural time and then analyzed in order to 
deduce the quantities ki, S and S- as well as the expo- 
nent aoFA from the slope of the log- log plot of the DFA 
analysis in natural time. We then examine whether the 
latter slope has a value close to unity or the conditions 
Ki « 0.070 and S, S- < Su are obeyed. In other words, 
the consequences caused by an undesirable severe data 
loss can be markedly reduced upon taking advantage of 
the DFA and natural time analysis. 
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FIG. 5: (color online) Electric field recordings at PAT in normalized units, i.e., by subtracting the mean value fi and dividing by 
the standard deviation a. The following SES activities are depicted in red: (a) on August 9, 2010, (b) on August 10, 2010. The 
dichotomous representations (blue) lead to the following values: ki = 0.072(4), S = 0.068(6), S- = 0.088(5) and ki — 0.064(5), 
S = 0.088(8), S- = 0.062(3), for (a) and (b) respectively. 




FIG. 6: The probability Prob(Ki) vs ki upon considering the seismicity (Mthre3=3.1) after August 10, 2010 until 15:42 UT on 
September 29, 2010 within the area A^ll.fSil.i. 




FIG. 8: Continuation of the study of the analysis of seismicity in natural time as described in the Note added on January 27, 
2011. 
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FIG. 9: The probability Prob(/ti) vs Ki upon considering the seismicity (Mt^T.es=l-9) after January 22, 2011 until 17:10 UT 
on March 7, 2011 within the area Aflfisll.'g. 




20°E 22°E 24°E 26°E 28°E 



FIG. 10: The area within which the seismicity was considered in natural time analysis after the initiation of the SES activity 
at PIR lasting from 23:00 UT on May 25 to around 02:00 on May 26, 2011. 
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FIG. 11: The probability Prob(«:i) vs ki upon considering the seismicity (Mthrcs=2.7) after the initiation of the SES activity 
at PIR (lasting from 23:00 UT on May 25 to around 02:00 on May 26, 2011) until 07:11 UT on June 24, 2011 within the area 
depicted in Fig. [TO] A maximum of Prob(/s:i) at ki ~ 0.070 is also observed for Mthres—'i.S upon the occurrence of the ML=4.2 
event at 06:53 UT on June 24, 2011. 
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*® Note added on April 16, 2010: Actually, as proposed 
in Ref. [82], three days later, i.e., on March 9, 2010, a 
Ms(ATH)=5.6 earthquake occurred at 38.87N 23.65E in- 
side the expected area. In addition, upon repeating the 
natural time analysis of seismicity*^ in the area N(37.65- 
39.0) E(22.2-24.1) (KER) after February 25, 2010, we again 
find an approach to criticality in the same fashion as de- 
scribed in Ref. [82]. The condition for criticality was found 
to hold only for magnitude thresholds Mthms = 2.8 and 
Mthres = 2.9 and hence in order to assure magnitude 
threshold invariance further study is in progress. 
Note added on May 29, 2010: The aforementioned study in 
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the area N(37.65-39.0) E(22. 2-24.1) finally revealed magni- 
tude tlireshold invariance of the cnticality condition upon 
the occurrence of the ML=3.0 event at 14:32 UT on May 
27, 2010, with an epicenter at 38. 5N 23.4E. This can be 
seen when plotting the probability Prob(Ki) versus ki for 
Mthres = 2.8, 2.9 and 3.0 in accordance with the procedure 
described in Ref. [85] identifying the occurrence time of an 
impending mainshock. 

Note added on August 23, 2010: The previous Note was 
followed on July 16, 2010 by Ms(ATH)=5.6 earthquake at 
39. 3N 24. OE, i.e., very close to the northeastern edge of the 
predicted area. Two other SES activities were recorded at 
Patras station (PAT) at 09:07 UT on August 9, 2010 and 
at 06:24 UT on August 10, 2010 depicted in Figs. 5(a) and 
5(b), respectively. They resulted in an analysis, in natural 
time, of the seismicity in the area N(37.6-39.0) E(20.4-23.2) 
similar to that in Ref. [82]. 

Note added on October 1, 2010: This is just a continu- 
ation of the analysis started, as mentioned in the previ- 
ous Note, in the area N(37.6-39.0) E(20.4-23.2) which is 
now extended to N(37.5-38.8) E(19.8-23.3). At 03:53 UT 
on September 2, 2010 a Ms(ATH)=4.8 earthquake actu- 
ally occurred with an epicenter at 38.22N 23.17E, i.e., in- 
side the area studied. In addition, upon the occurrence of 
a Ml =3.1 event at 15:42 UT on September 29, 2010 with 
an epicenter at 37.56N 20.00E, when plotting the proba- 
bility Prob(Ki) versus ki, for Mt(ir-es=3.1, we find a sharp 
maximum at ki ~ 0.070 depicted in Fig[6l which probably 
indicates that the critical point has been approached. To 
ensure that this is a true criticaltty condition, the extent 
to which a magnitude threshold invariance holds for this 
result is currently investigated. 

Note added on October 28, 2010: Actually, almost one 
week after the previous Note, i.e., at 19:04 UT on Octo- 
ber 9, 2010, a Ms(ATH)=4.8 earthquake occurred with an 
epicenter at 38.15N 22.72E lying inside the area studied. 
The continuation of the study related to the SES activi- 
ties recorded at PAT on August 9 & 10, 2010 (FigO in 
the same area N(37.5-38.8) E(19. 8-23.3), reveals the fol- 
lowing: Upon the occurrence of a ML=4.1 earthquake at 
10:32 UT on October 23, 2010, with an epicenter at 38.73N 
21.99E, a maximum in the plot of Prob(«:i) versus ki for 
Mthi es=3.3 is observed at ki = 0.070. The same was ob- 
served for Mtfir(!a=3.2 upon the occurrence of a ML=3.2 
earthquake at 00:38 UT on October 24, 2010, with an epi- 
center at 38.71N 21.99E as well as for Mthres=3.1 upon the 
occurrence of a ML=3.6 earthquake at 04:04 UT on Octo- 
ber 28, 2010, with an epicenter at 38.36N 22.25E. These 
findings are strikingly reminiscent of the following behav- 
ior: For higher magnitude threshold, the description of 
the real situation approaching cnticality becomes less ac- 
curate due to "coarse graining" . The plot for the aforemen- 
tioned maximum observed on October 28, 2010 is depicted 
in FiglZl 

Note added on January 27, 2011: At ^ 01:00 UT on 



January 22, 2011 an electric disturbance started at KER 
which lasted for around 9 hours. This does not obey the 
Ay/L criterion so that to be classified as SES activity. 
On the other hand, when analyzing it in natural time we 
find the following values ki = 0.072(6), 5" = 0.075(6) and 
S- — 0.085(9) which conform with the properties of a true 
precursory signal. A violation however of the AV/L cri- 
terion may be due to the inhomogeneities in that area, 
and in order to resolve the possibility of an approach to 
criticality, we continued the investigation described in the 
previous Note as well as in the area A^ll.l Eii'.e,- We find 
that (for Mthres=3.2) upon the occurrence of a ML=3.2 
event at 09:31 UT on January 24, 2011, a maximum of 
Prob(Ki) versus ki is observed at ki = 0.070, see Fig[S] 
which possibly indicates the approach of the system to the 
critical point. 

Note added on March 8, 2011: Upon the occurrence of the 
ML=3.5 event at 17:10 UT on March 7, 2011, the Prob(Ki) 
versus ki exhibited a maximum at the value ki = 0.070 
for Uthres^l-l to 1.9 (see Fig. |9] for Mth,e«=1.9). If the 
seismic data for such low magnitude thresholds are com- 
plete (which is questionable) , the present behavior suggests 
that the critical point is approached (i.e., the event is ex- 
pected within a few days to around 1 week or so) in the 
area N^iElti- Since AV/L is of the order of 5mV/km, 
Ms(ATH) depending on the epicentral distance may reach 
Gifr ~ 100km. 

Note added on June 24, 2011: Since the seismic data 
in the previous Note for such low magnitude thresholds 
were not complete, as questioned, the time-window calcu- 
lation could not be accurate. An earthquake of magnitude 
Ms(ATII)=4.4 actually occurred later on AprU 22, 2011 
with an epicenter at 38.37°N 23.62°E, i.e., around 60-70km 
from KER. In addition, a SES activity was recorded at 
PIR lasting for around 3 hours, i.e., from 23:00 UT on 
May 25 to around 02:00 on May 26, 2011. Thus, a nat- 
ural time analysis of the subsequent seismicity started in 
the shaded area depicted in Fig llOl in a similar fashion as 
in arXiv:0802.3329v4 (publicized on May 29, 2008 before 
the 6.5 on June 8, 2008) following the same procedure de- 
scribed in the latter publication. This revealed that, upon 
the occurrence of a ML=2.9 earthquake at 07:11 UT on 
June 24, 2011 with an epicenter at 37.6°N 21.0°E, the 
Prob(«:i) versus exhibited a maximum at ki « 0.070 
(as shown in Fig. Ill[). thus pointing to the approach of the 
system to the critical point. 

Note added on January 26, 2012: In a recent manuscript 
submitted for publication on January 23, 2012, the 
recording of a strong electric disturbance at PIR (see 
Fig llOp on January 6, 2012 has been reported. A copy 
of its recording along with the current investigation on 
its tentative classification as SES activity can be found 
elsewhere (see the file arxiv_230112.pdf available from 
|http: //physlab.phys .uoa.gr/org/director.html ). 



